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Abstract 

Based on the method of collective variables we develop the statistical field theory for the study 
of a simple charge-asymmetric 1 : z primitive model (SPM). It is shown that the well-known 
approximations for the free energy, in particular DHLL and ORPA, can be obtained within the 
framework of this theory. In order to study the gas-liquid critical point of SPM we propose the 
method for the calculation of chemical potential conjugate to the total number density which 
allows us to take into account the higher order fluctuation effects. As a result, the gas-liquid phase 
diagrams are calculated for z = 2 — 4. The results demonstrate the qualitative agreement with 
MC simulation data: critical temperature decreases when z increases and critical density increases 
rapidly with z. 



I. INTRODUCTION 



Besides being of fundamental interest, ionic systems including electrolyte solutions, 
molten salts and ionic liquids deserve great attention from practical point of view. For 
example, new ionic liquids with very low vapor pressure may find applications in environ- 
mentally clean industry. Over the last ten years the phase and critical behavior of ionic 
fluids has been a subject of intense research. For reviews of experimental and theoretical 
situation see Refs. and the references cited therein. Basic properties 

of ionic systems can be described by the primitive model (PM) consisting of a mixture of 
m species of charged hard spheres. A two-component mixture of positive go an d negative 
— zqo charges having all the same diameter a is a simple charge-asymmetric PM (SPM). In 
the case z = 1 we have a well-known restricted primitive model (RPM). Early studies ^] 
established that the RPM has a gas-liquid (GL)-like critical point. A reasonable theoretical 
description of the GL critical point of the RPM was accomplished at a mean-field (MF) 
level using integral equation methods and Debye-Hiickel (DH) theory |lOj. Over the 

last decade numerous simulation studies have been devoted to the location of GL phase 



transition of this model 



and the most reliable current 



estimations turn out to be near T* = 0.049, p* = 0.08 when the temperature T* and the 
density p* are in standard dimensionless units. Due to controversial experimental findings, 
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the critical behavior of the RPM has also been under active debates 
strong evidence for an Ising universal class has been found by recent simulations 
and theoretical |28 |29| studies. 

In spite of significant progress in this field, ionic systems are far from being completely 
understood. The investigation of more complex models is very important in understanding 
the nature of critical and phase behavior of real ionic fluids which demonstrate both the 
charge and size asymmetry as well as other complexities such as short-range attraction. The 
studies of the effects of charge-asymmetrv on the phase diagram have been recently started 
usmg both th e c 0mpute slmul a tl o ns fly md thecal — WW- 
Monte Carlo (MC) simulations show that SPM exhibits a GL-like phase transition, specifi- 
cally T*(z) decreases with the increase of z while p* c {z) rapidly increases. Unfortunately, the 
results obtained within the framework of both the original DH theory and the mean spher- 
ical approximation (MSA) are independent of z. The results for the z-dependence of the 
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critical temperature, obtained within the framework of the field-theoretical analysis 
contradict the simulations displaying the increase of temperature with z. On the other hand, 
the recent results found based on the DH theory incorporated Bjerrrum ion pairs and their 
solvation in the residual ionic fluid j^H^I (theories DHBjCI and DHBjCIHC) demonstrate 
the decrease of T* with charge asymmetry for z = 1 — 3 which agrees with MC data. This 
implies "... that recognizing ionic association is inescapable for a successful theory" 

The motivation of this paper is to check the above statement by means of the theory that 

n 

exploites the method of collective variables (CVs) 36]. Our aim is to answer the question: is 
the theory, which does not include the effects of association explicitely in the Hamiltonian, 
capable of correctly describing, at least on the qualitative level, the GL-like phase diagram 
of SPM? To this end, we develop a theoretical approach which is based on the functional 
representation of configurational Boltzmann factor in terms of CVs. For a general case of a 
two-component model which includes both the short- and long-range interactions we derive 
an exact expression for the functional integral of the grand partition function (GPF). The 
CV action obtained depends upon the two sets of CVs: {pk} (and conjugate {o>k}) and {Qk} 
(and conjugate {7k}) describing the total number density and charge density fluctuations, 
respectively. We start with the saddle point (or MF) solutions p (and uj) and Q (and 7), 
and then we expand the CV action functionally around the MF solutions. As a result, 
we obtain an infinite expansion in terms of CVs <5pk and 5Q^ (or Su^ and 57k). For SPM, 
based on the GPF in the Gaussian approximation, we propose a method for the calculation of 
chemical potential conjugate to the total number density. The method allows us to take into 
account the higher order fluctuation effects and consists in solving the equations for chemical 
potentials by means of successive approximations that correspond to the contributions of 
the higher-order correlations. Its initial idea to some extent resembles the idea sketched out 
by Hubbard in [37 1. 

It should be noted that the CVs were first introduced in the 1950s and then 

on their basis the so-called CVs method was developed 0, 3|. Nearly at that time, 



other functional approaches based on the Stratonovich-Hubbard transformation 



were originated. As was shown recently 



HQ 



42j both groups of theories are in close relation. 



The layout of the paper is as follows. In Section 2, starting from the Hamiltonian of an 
asymmetric two-component model with long- and short-range interactions, we derive the 
exact functional representation of the GPF. We obtain an explicit expression for free energy 
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of SPM in the random phase approximation (RPA) that leads to the well-known results, 
i.e. the DH limiting law (DHLL) for the point charge particles and the free energy in the 
optimized random phase approximation (ORPA) (or MSA) for SPM. In Section 3 we study 
the GL critical point of SPM taking into account the correlation effects of higher order. The 
results obtained demonstrate qualitative agreement with the MC simulation data for both 
the critical temperature and the critical density. We conclude in Section 4. 

II. FUNCTIONAL REPRESENTATION OF THE GRAND PARTITION FUNC- 
TION OF A TWO-COMPONENT IONIC MODEL IN THE METHOD OF CVS 

Let us consider a general case of a classical two-component system consisting of N par- 
ticles among which there exist N± particles of species 1 and N2 particles of species 2. The 
pair interaction potential is assumed to be of the following form: 

M^W + ^W + W. (i) 

where <f>ap( r ) is the interaction potential between the two additive hard spheres of diameters 
a aa and a pp. We call the two-component hard sphere system a reference system (RS). 
Thermodynamic and structural properties of RS are assumed to be known. (f)^p(r) is the 
Coulomb potential: <fiap( r ) — QaQp<p C '(r) , where c (r) = l/(Dr), D is the dielectric constant, 
hereafter we put D — 1. The solution is made of both positive and negative ions so that the 
electroneutrality is satisfied, Y^ a =i QaC a = 0, and c a is the concentration of the species a, 
c a = N a /N. The ions of the species a — 1 are characterized by their hard sphere diameter 
(Tn an d their electrostatic charge +q and those of species a — 2, characterized by diameter 
(T22, bear opposite charge — zq (q is elementary charge and z is the parameter of charge 
asymmetry). 0f^(r) is the potential of the short-range interaction: <fr%p(r) = 4>ap{ r )+4>tp( r )^ 
where (fiap( r ) * s use d to mimic the soft core asymmetric repulsive interaction, 0fjg(V) is 
assumed to have a Fourier transform; 4>ap{ r ) describes a van der Waals-like attraction. 

We consider the grand partition function (GPF) of the system which can be written as 
follows: 



=2^ 2^ [[ — J ( dr ) ex P 



aP ij 



(2) 



JVi>0 JV 2 >0a=l,2 

Here the following notations are used: v a is the dimensionless chemical potential, v a = 
Pfj, a — 3 In A, ji a is the chemical potential of the ath species, j3 is the reciprocal temperature, 



A -1 = {2nm a f3~ l /h 2 ) 1 / 2 is the inverse de Broglie thermal wavelength; (dr) is the element 
of configurational space of the particles. 
Let us introduce operators /5k and Qk 



which are combinations of the Fourier transforms of the microscopic number density of the 
species a: p^ a = ^ exp(— ikr"). In this case a part of the Boltzmann factor in (J2J) which 
does not include the RS interaction can be presented as follows: 



cxp 



-§££(^«)-0*)) 



a/3 i,j 



cxp 



■5E(* 



NNPkf>-k 



(3) 



where 



i 

i 



z 2 $fi(k) + 2z$ s 12 R (k) + $ s 22 R (k) 
®f?(k) - 2$ff(Jfe) + <£|f (&)] + $ c (k) 
z*f*{k) + (1 - z)^{k) - $$*(kj 



(4) 



and we use the notations = y4%B'{k) with 0^g"(fc) being a Fourier transform of the 

corresponding interaction potential. 

In order to introduce the collective variables (CVs) we use the identity 



exp 



l . r i - 

-2j$W7- k = (drj)5r[ri-fj\exp --zz2®Wn 

l, ** l. 



(5) 



where — 77] denotes the functional delta function 



fcfo - »)] = / ( dA ) ex P 



i£A k (r7 



?7kJ 



= "^k — i^k (^k = Pk, <5k) is the collective variable and 

(dr?) = dr/o J] dt£di£, (dA) = dA ]J dA k dA k . 

k^O k^O 

The indices c and s denote real and imaginary parts of r/k (Ak), respectively, the product 
over k is performed in the upper semi-space. 
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Taking into account we can rewrite (J2J) 

S M = J (dp)(dQ)(d^)(d7)exp(-ft[z/ a ,p,Q,w,7]), (6) 
where the CV action 7i is as follows: 

H[v a , p, Q, co, 7] = - ^[$ivjv(A;)pkP-k + ®QQ(k)Q k Q- k + 2§ NQ (k) 

k 

xpkQ-k] - iJ^kPk + 7kQk) - \nE HS [u a ; -iu, -iq a j]. (7) 

k 

In ((7|) CVs pk and Qk describe fluctuations of the total number density and charge density, 
respectively. Ens[^a', —iui, —iq a l\ is the GPF of a two- component system of hard spheres 
with the renormalized chemical potential v a in the presence of the local field ip a {fi) 

exp(u a N a 



s»[-] = ee n- 

7Vi>0 7V 2 >0a=l,2 

-EE^( r * 



(dr) exp 



kHSi 



9 EE<w r ^ 



a/3 ij 



where 



£E *£(*) + # E*°(*)' 



(8) 



(9) 



Vfefa) = iu(ri) + \q a l{ri). (10) 
a. Mean-field approximation. The MF approximation of functional (jUJ) is defined by 



SmfM = exp(-W[f a ,p,Q,a;,7]), 
where p, Q, u and 7 are the solutions of the saddle point equations: 

p = {N[u a ; -iu>, -iqa^HS, = 0, 
u = -ip4> NN (0), 7 = -ip4> NQ (0). 

Substituting (fT2"j) in (fTT|) we obtain 



(11) 



(12) 



^mfM = exp 



L ^P 2 0w7v(O) 



z/ a ; -iu;,-ig Q 7]. 



The Legendre transform of In S m p gives the free energy in the MF approximation 

A 

2V 

2' 



2^ 



(13) 



It is worth noting that for </>£g (r) = we arrive at the free energy of SPM in a zero- loop 



approximation 

b. Beyond the MF approximation. First we present CVs p k and Q k (w k and 7k) as 

Pk = pS k + 8f>k, Qk = Q$k + SQ k , 
u k = td5 k + Su; k , 7k = 74 + <$7k, 



where the quantaties with a bar are given by (jl2j) and 5 k is the Kronecker symbol. 

Then we write In E H g[p a ; — iui, — iq a j] (see (jH|)- (fTU|) ) in the form of the cumulant expansion 



\nE HS [...} = E W$ n) fa,...,k n )6y kl ...&r kin 

n>0 i n >0 ki,...,k n 

Su *i n+1 ■ ■ ■ Suj kJk 1 +.„+k n , (14) 
where DJtf£ n \ki, . . . , k n ) is the nth cumulant defined by 

^['"'(fcb ■ ■ ■ ,k n ) = T7T 7T? KT~~ ■ ^7 |5 7k =0,<W k =0- (15) 

d5j kl . . . d5j kin dSuj kin+i . . . d5u kn 

The expressions for the cumulants (for n < 4) are given in Appendix A. 
We can integrate in (jUJ) over 5u k and #7k and obtain for E[u a } (see j^j) 



-- EmfWoI^' J (d5p)(d5Q) exp | - ^ ^ [L NN (u a ; k)5p k 5p^ k 
J ' k 

+2L NQ (u a ; k)Sp k SQ^ k + L QQ (v a ; k)5Q k 5Q_ k ] 

+EE^ n) (^^^)}- ( 16 ) 

n>3 i n >0 

It is worth noting that the above expression is of the same form as that obtained within the 
framework of the mesoscopic field theory j^l^. The main difference is that H[z/ a ] in ()16j) is 
a function of full chemical potentials, rather than just a function of their mean field parts. 

A. Gaussian approximation. A two- component primitive model (PM) 

In the Gaussian approximation, which corresponds to taking into account in (j!6|) only 
terms with n < 2 (Tin =0), we have Lab^u] k) = Cab^u] k) (A, B = N,Q), where 



C NN = $ NN (k) + l/lf(i/ a ; k), C QQ = $ QQ (k) + l/mP(u, 



a J i 



C NQ = $ NQ (k). (17) 



Cjvjv(^), Cqq(^) an d CjvQ(fc) are the Fourier transforms of the density- density, charge-charge 
and charge-density direct correlation functions, respectively. It is worth noting that DJt^ {v a ) 
does not depend on k (see Appendix A). After integrating in (jTfij) over <5p k and 5Q^ taking 
into account (|17j) we arrive at GPF in the Gaussian approximation. 

Let us consider a two-component simple primitive model (SPM) j3] consisting of charged 
hard spheres of the same diameter {an = o~22 = cr) with $£g (&) = which differ by their 
respective charges (z 7^ 1). We have for SPM (see 

®NN(k) = 0, § NQ (k) = 0, $ QQ (A;) = i> c (A;). (18) 

It should be noted here that the Hamiltonian of SPM, as for RPM 2^| , does not include 
a direct pair interaction of number density fluctuations. Integration over <5pk and Slv-^ in 
(JHJ) is trivial in this case and leads to the KSSHE action obtained in [143] by performing the 
Hubbard-Stratonovich transformation. Starting from this expression the free energy of SPM 
in a two-loop order approximation was derived by J.-M. Caillol As was shown jss] ]. the 
z dependence of critical temperature calculated within the framework of this approximation 
disagrees with the results obtained by MC simulations. Below we propose an alternate way 
of taking into account the fluctuation effects near the GL critical point. 

Expressions flT7|) . under conditions (|TH|). are reduced to 

C NN = 1/Wl { °\v a ; k), C QQ = $ c (k) + lrf(^), C NQ = (19) 
and the logarithm of GPF in the Gaussian approximation is as follows: 

\nE G (u a ) = In ~hs(^) -\j2 ln [ 1 + ®° \ k )^ (*«)] ■ ( 20 ) 

k 

After the Legendre transform of lnH<3(z/ a ) we obtain the well-known expression (see e.g. 
4^ |) for the free energy in RPA 

PIrpa = PfMF + ^7 J>(1 + K 2 D <f> c (k)), (21) 

k 

where fifuF has the form (fT3")l under conditions (fTSj) . k 2 d = 4irpj3q 2 z is the squared Debye 
number. It is worth noting that a use of momentum cutoff |1ca| = 2-K/a in (J2"T|) leads to 



the same expression for the PJdh as in 



3- 



For point charge particles, ()21|) yields the free 



energy in the DHLL approximation (3fDHLL = Pfid — j^- Using the optimized regularization 
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of the Coulomb potential inside the hard core [47( we arrive at the free energy in ORPA (or 
MSA). 

As is seen, P/rpa does not explicitely depend on the charge asymmetry factor z. The 
same is true for MSA and the DH theories. The detailed analysis of the GL phase equilibrium 
using (|21|) and different regularizations of the Coulomb potential inside the hard core was 
fulfilled in j^J. As was shown 3^, the ^-dependent free energy can be found only in 
the higher order approximations. 

Below we study the GL phase diagram of SPMs when fluctuation effects of the order 
higher than the second order are taken into account. 

III. GAS-LIQUID CRITICAL POINT OF THE PRIMITIVE MODELS (PMS) OF 
IONIC FLUIDS 

Let us consider equation (f2~U|) . Introducing v^ and vq 

" N = ih" i+ TTz^ * q = wvz)^-^ (22) 

we can present (J20J) as follows: 

ln~ G K) = In ~ hs (^q) + * C fa > *q)] > ( 23 ) 

k 

where new chemical potentials v^ and uq (see (j22|0 are conjugate to the total number 
density and charge density, respectively. Since near the GL critical point the fluctuations of 
the number density play a crucial role, v^ is of special interest in this study. 
We present v N and vq as 

V N = V% + \l±V N , Vq = V Q + XAv Q , 

with z/Jy and v Q being the MF values of v^ and vq, respectively and Av^ and Avq being 
the solutions of the equations 

d In E G (v N , v Q ) 



dAv N 
d In E G (v N , v Q ) 



dAv Q 

Expanding ()23|) in powers of Av^ and Avq we obtain 



X(N), (24) 
0. (25) 



n A/f (in) (z/ v°) 

In E G (u N , v Q )=J2I2 <% H n ^N in i (26) 



n>0 i„>0 



where 



KA{jn) (l P y Q) - d n \nE G (is N ,is Q ) 



The expressions for the coefficients A4n are given in Appendix B. 

We solve equations (j23)-(I2ni) approximately taking into account and keeping terms 
of a certain order in parameter A. The procedure is as follows. First, we calculate Auq 
from (J2I3) in the approximation which corresponds to a certain order of A e.g., order s. 
Then, we substitute this Auq into equation (}2"4"j) in order to find Az/jv in the approximation 
corresponding to A s+1 . In ()24)1 we take into account only the linear terms with respect to 
Ai/jv keeping terms with all powers of Aug within a given approximation in terms of A. 

As is readily seen, the first nontrivial solution for Az//y is obtained in the approximation 
of the first order of A. It is the result of substitution in (|2*4"j) of the solution Apq = 0. As a 
result, we have 

where 

S(k) = ® C(k) = 1 (28) 

(f2Tj) - (j2*Hj) can be rewritten as (see Appendix C) 

Az/jv = > D ~ 1 29 

which corresponds to RPA. As was mentioned above, Av^ given by (}2*9~|) does not depend on 
charge asymmetry factor z. In order to obtain the ^-dependent expression for the chemical 
potential related to the number density fluctuations we should consider the next approxi- 
mation in A for Av^. 

In order to obtain Az/jy in the approximation corresponding to A 2 we substitute in (I24|) 
the solution Auq as follows (see Appendix C): 

= -|jj&£sr(*) = £«*(*). (30) 

ZJJl 2 k k 
which is found from (|25jl in the first approximation of A. Taking into account only a linear 
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term with respect to Av^ we get 

1 

/\u N = 



m 2 



(31) 



(0) 

1 (9\ . r, 1 



2 ' 6 2 

k k 



Let us consider (|31|) in detail. The correlation effects of the order higher than the second 
order enter the equation through the cumulants 9Jln™' ) for n > 3 and i n ^ 0. The appearence 
of these cumulants reflects the fact that the terms proportional to wy 2 , 7 s and or/ 3 are 
taking into account in the cumulant expansion ()14|) (n < 4). Recall that and 7k are 
conjugate to the CVs pk and Qk describing the total number density and charge density 
fluctuations, respectively. This means that in order to determine Az//v we take into account 
in (|14j) . besides the terms of the second order, the contribution corresponding to the pure 
charge density fluctuations for n = 3 (JXH 3 ) and the contributions corresponding to the 
correlations between charge density and total number density fluctuations for n < 4 (9Jt^ 
and OJQ ) which are linear in u^. Therefore, the analysis of (}3*T]) establishes a link between 
the approximations considered above (in terms of A) and the approximations formulated in 
terms of CVs. 

Another important issue to be discussed is the limitting case z = 1 that cooresponds 
to RPM. For z = 1 only the first term survives in (|31j) . Furthermore, in this case the 
conditions DJln^ = and Auq = hold simultaneously (they are equivalent for the RPM) 
due to the model symmetry (see formulas in appendixes A and C). As a result, we arrive at 



33, 



481 ] . On the other hand, when 



the expression for the chemical potential of RPM in RPA 
we set solely Az/q = in (|3*T|) we obtain (|2T|) which also corresponds to RPA and is valid for 
z 7^ 1. The latter reflects the fact that the correlations between fluctuations of the charge 
density and the total number density are not taken into account at the RPA level. Therefore, 
for Az/q = all the terms, except the first one, in (}3Tj) become equal to zero and we get the 
expressions for the chemical potential conjugate to the total number density in the same 
approximation for the both models. The equation (}3T|) is obtained in the approximation 
when only the liner terms with respect to uo in the cumulant expansion (fT4*|) are taken into 
account. However, for z = 1 these contributions are equal to zero. This means that the 
approximation considered in this paper for SPM does not have an analogy for RPM. In order 
to include the fluctuations in the simplest model (RPM) the higher order terms should be 
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in 3 



where a good agreement for the critical 



taken into account. This issue was considered 
temperature was obtained. 

Based on (|31|) the coexistence curves for SPM for differnt values of z can be calculated. 
It is worth noting that the regularization of the potential inside the hard core is 

arbitrary to some extent. For example, different regularizations for the Coulomb potential 
were considered in jf], . Within the framework of the Gaussian approximation of GPF 
;he best estimation for the critical temperature is achieved for the optimized regularization 



3 



47| that leads to the ORPA (MSA). However^ 



lis approximation does not work properly in 
. In this study we use the Weeks- Chandler- 



the higher orders of the perturbation theory 

Andersen (WCA) regularization j^j]. As was shown 46j . this choice of cf) c {r) for r < o 
produces rapid convergence of the series of the perturbation theory for the free energy. For 
the WCA, the Fourier transform of <j)°(r) is of the form (f> (x) = sin(a;)/x 3 with x = ka. As 
a result, (|3"Tj) can be written as (see Appendix C) 



7i 



2zti 



3zir 



In 



x 2 sin xdx 



(32) 



(33) 



jo J " 3 

with k 

introduced for the temperature and the density: T" = , p" = pa 

Finally, we get the expression for the full chemical potential conjugate to the total 
number density 

z 



o x d + k* smx 

K D a being the reduced Debye number. Hereafter the standard notations are 

k B Tcr * 3 



vn — 3 In A/cr 



. 77(8 - 9?7 + 3t7 3 ; 
mp H h 



0--VY 



l + z 



In z — ln(l + z) 



2T* 



(34) 



where Az/jy is given in ()32j) - (f3"3|) and r\ = Tip* /Q. In ()34jl the Carnahan-Starling approxima- 
tion for the hard sphere system is used. 

Fig. 3 shows the coexistence curves for 2 = 2 — 4 calculated based on the isotherms of 
chemical potential (pMj) supplemented with the Maxwell construction. The estimations for 
the critical point are as follows: 

z = 2 : T* = 0.12310, p* = 0.00946 
z = 3 : T* = 0.11313, p* = 0.02740 
z = 4 : T* = 0.10030, p* = 0.04501. 
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FIG. 1: Coexistence curves of the (1 : z) charge-asymmetric ionic model. 

These results demonstrate the qualitative agreement with the MC data: the critical temper- 
ature decreases when z increases and the critical density increases rapidly with z. Moreover, 
a comparison of the coexistence curve forms for z = 2 and z = 3 (results for the coexistence 
curve for z = 4 are not available by now) with the DHBjCINC theory indicates their simi- 
larity. It is also worth noting that the above data for the critical temperature lie about in 
the same region as those obtained in j3| within the framework of the two-loop expansion. 




IV. CONCLUSIONS 



In this paper we have studied the GL phase behaviour of a charge asymmetric primitive 
ionic model. For this purpose we have derived the exact expression for the functional of GPF 
of a two-component asymmetric ionic model which includes both the short- and long-range 
interactions among charged hard spheres. We have shown that the well-known approxi- 
mations for the free energy, in particular DHLL and ORPA, can be reproduced within the 
framework of this approach. On the other hand, the GPF functional can be reduced to the 
form found in the KSSHE theory This means that the field-theoretical analysis of the 
expression for GPF given by ©-© has to lead in the two-loop approximation to the same 
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z-dependence for the critical temperatures as in |33J which does not agree with the MC 
simulations. Here, we have proposed an alternative method for the study of the GL phase 
equilibria in SPM. It consists in the calculation of the chemical potential conjugate to 
the total number density by means of successive approximations. 

In conclusion, we have obtained the expression for the chemical potential in which 
the effects of indirect correlations between the number density fluctuations are taken into 
consideration via a charge subsystem. Based on this expression supplemented with the 
Maxwell constraction the coexistence curves for ,2 = 2 — 4 have been calculated. The z- 
dependences obtained for both the critical temperature and the critical density qualitatively 
agree with MC simulations. The results demonstrate that the terms responsible for the 
charged clustering need not be included into the Hamiltonian explicitely in order to describe, 
at least qualitatively, the GL phase diagram of SPM. Instead, the interaction between charge 
and density fluctuations should be properly taken into account. To achieve a quantitative 
agreement, the higher order terms should be taken into consideration in the CV action. This 
task will be considered elsewhere. 
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Appendix A: Recurrence formulas for the cumulants in the Fourier space represen- 
tation 



m^ih,^,...,^) = G n {h,k 2 ,...,k n ) 
rni^(k h k 2 ,...,k n ) = o 

9ft%\k 1 ,k 2 ,...,k n ) = ^ qlc a G n „ 1 (k 1 ,k 2 , ...,1k, 



■n— 



i + K\) 



(35) 
(36) 
(37) 



a 



m^ih, k 2 , . . . , k n ) = J2<ll c aGn-2{ki,k 2 ,... 



k n _ 2 + k n _i + k„|) 



(38) 



a 
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where GJb,h,...,K> * the Fourier transform of the particle truncated correlation 
function |4J] of a one-component hard sphere system. 



Appendix B: Explicit expressions for A4n 





= x^ + ^^k), 






M? 


k 

k 






MP 


= x^ + ^m^ 

k 


A 6 

+ y 


k 




k 




k 






A 6 

+ y 





(0) 



(2) 



k k 

+X*[mf\vlv° Q )] 3 J2~9 3 (k), 

k 

£m?> e^) + t K< } + e*w 



+A 9 [^ 2) ]V ) E^( fc )' 

k 

A 3 ^ 2) + ^ 4) E^) + y 

k 



xanfK,^)jE^) + A £ 

k 

E^)> 



n 2 



X 



(3) 



A 3 ajtf + £m?> E^) + ^ 3) *ti 4) E^ 2 W 



+A fe 



*4 3) 1 E^ 
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M 



(0) 



M 



(i) 



M 



(2) 



if 



(2) 



A 4 ^ + ^ 2) ^^) + ^(3[^ 

k 

x 9 2 (k) + 6A 10 [£E4 2) ] 'mtf ^ ^ 3 (fc) + 3A 12 071 

k k 



2 7 2 

k 

< ^ 2 (ifc) + 3A 10 ^fanfmif + 



971 



(2) 



+3A 12 W 2) 



+mtl 2) 27ti 4) + 2 [sjtf] 2 + 2971? [jmf ] 2 ) £ 5 a (Jfe) + 2A 

k 



x ( 2Wlfwif9Jlf + 



97t 



X 



A J 



971 



(3) 



m f) + £ g(k) + ^ (sanfmtf + mf m. 



A _! 

2 V 



If 



+3Tlfmf) ^~g 2 {k) + 3A 10 (97lf!mf 971, 



art: 



(3) 



n 2 



97t 



(3) 



^^(fc) + 3A 12 97tr ) 97tf ^^ 4 (fc). 



In the above formulas 97Tn 



(in) 



Appendix C: Some explicit relations obtained for a 1 : z asymmetric PM 



Let us consider the expressions ()35J ) - (j3*Bj) . For 1 : z asymmetric PM we have 
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Using the above relations and (|35 |) -(|38 |) we can obtain the following explicit relations for 
SPM 

gffi _ 2 V _ 2 (i - zf mfmf _ 2 (i-*) 1 





2 




= & 
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